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I. INTRODUCTION 


A.s SObUBCEIV: 

This thesis is an investigation of synthetic aperture 
radar (SAR) imagery as a tool for hydrographic surveying. 
The objective is to relate radar backscatter from surface 
wave phenomena to bathymetry and ocean bottom features. A 
descriptive interpretation of macro-scale features, 
(including internal wave signatures which are visible to the 
eye), and pseudo-imaging of shoals, is analyzed. Also, the 
Optical Fourier Transforms (OFT's) are used to estimate 
surface gravity wave refraction from SAR images of selected 
areas in Shelikof Strait, Alaska; linear wave refraction 
estimates are then used, in turn, to estimate water depths 


in shoaling areas. 


B. BACKGROUND 

Traditional waterways were established for ships with 
drafts of 35 feet and less, and are inadequate for modern 
containerized shipping methods and tankers (Goldsteen, 
1982). Therefore, new shipping lanes must be found to allow 
the new, deeper-draft vessels to navigate waters which were 
previously unused. However, conventional hydrographic 


resources are inadequate to provide the accurate nautical 
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Charts now required by both civilian and military users of 
these newer, larger ships (Hammack, 1977). 

Only 16 percent of the world's oceans has adequate 
Survey data to accurately represent the ocean bottom 
topography, and the time required to survey the remaining 84 
percent using conventional survey methods is enormous. For 
example,-a British study concluded that 284 ship years of 
work would be required to adequately survey Great Britain's 
territorial waters.  (Hammack, 1977) 

There are also severe financial constraints associated 
with acquiring detailed survey data. A well-equipped 
coastal survey vessel costs approximately $20 million and 
has yearly operating expenses of $1.4-3.0 million 
(Goodfellow, 1982). The cost of an average survey approaches 
/$3500 per sque mile (Collins, 1978). The Shelikof Strait 
areas examined in this study cover approximately 7,000 
square miles. Thus, it would cost $24.5 million to 
adequately survey this area. 

Because survey resources are limited (Hammack, 1977), a 
remote sensing instrument, such as SAR, could be used as a 
planning tool for directing the application of those 
resources at a minimal cost. The imagery used in this study 
cost $2,500 for processing; and the cost of collecting the 
raw data should be negligible once the satellite is in 


orbit. 
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A synthetic aperture radar measures a signal 
backscattered by the Bragg scattering cross section of 30 cm 
ultragravity waves on the ocean surface. Bathymetric 
features have been shown to influence this cross section 
indirectly, with the effect of revealing the shapes of 
shoals in distinct patterns of brightness modulation in the 
‘SAR image (Plant, et al., 1978, Evans and Shemdin, 1980, 
Alpers, et al., 1981, Brown, et al., 1976, McLeish, et al., 
1980 and Kasischke, et al., 1983). Perhaps the most 
dramatic example is the almost one-to-one relationship 
between the SAR-imaged ultragravity waves and the bottom 
sand waves in Nantucket Shoals (Kasischke, et al, 1980 and 
Valenzuela, et al., 1983). This example will be reviewed in 
Chapter II of this thesis. 

Other examples of oceanic processes that have distinct 
SAR signatures are: (1) wave refraction, (Shuchman, et al., 
1979), (2) non-linear gravity wave interaction, (Kasischke, 
et al., 1980, (3) upwelling, (Kastscehke, ect al.) eae ans 
(4) internal E manifestations (Shuchman and Kasischke, 
1979). Wave refraction occurs as a surface gravity wave 
propagates shoreward over a shoaling bottom; both the wave 
direction and wavelength change as a linear function of 
water depth. If a surface gravity wave encounters an abrupt 
change in depth, such as a coral reef or submerged rock, the 
ensuing wave transition strains the ultragravity wave field. 


Colder, deep water rising to the surface (upwelling) 
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modulates the ultragravity wave structure in two ways: 1) 
the convergence of temperature and biological influences 
vary surface tension as a restoring force, and 2) the 
planetary boundary layer is more stable over colder water 
than over warmer water, and hence, wind stress is lower, 
showing up as darker areas in SAR imagery. When a current 
flows over a topographic bottom feature, internal waves can 
be generated which modulate the surface ultragravity wave 
field. 

The evidence is clear that SAR signatures indicate the 
presence of various bottom configurations, and that SAR can 
be used to locate at least some areas of potential danger to 
ocean and coastal navigation. SAR signatures of bathymetric 
phenomena may be catologued for use by an image analyst, and 
wave refraction estimation techniques may be used to 
estimate depths from SAR imagery. A hydrographic survey 
planner can use this type of information to optimize 
allocation of survey resources, scheduling high priority 
areas where possible hazards are indicated for immediate 
surveying, and leaving other areas for surveying at a later 


date. 


Gye PROACH 
Shelikof Strait, Alaska is typical of areas of 
consideration for new shipping routes, and extensive 


resources would be required to adequately survey the area. 
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1549 


Figure 1.1. SEASAT Coverage ninel. Sa a Ss a 
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Therefore, SAR images were obtained for analysis and 
pumuvmetric maps of Shelikof Strart" were obtained to usé as 
adl reference for the analysis. 

The SAR images and subareas were first examined to 
identify obvious anomalies or patterns in the radar cross 
section modulation which might be associated with 
bathymetry. These signatures were then compared to colocated 
bathymetric features as shown on bathymetric maps of the 
area. The bathymetric maps were also reviewed to identify 
areas where gravity wave refraction was likely to occur. In 
those selected areas, Optical Fourier Transforms (OFT) were 
acquired and analyzed to infer surface gravity wave 
refraction (change in wavelength and direction) for 
comparison with changes predicted using linear wave 
refraction and charted water depths. 

Shelikof Strait was chosen as the study area due to the 
availability of both bathymetric maps based on recent 
surveys and SAR imagery from three SEASAT orbits (Fig. 1.1). 
Also, this area is representative of the types of waterways 
which are candidates for new shipping routes. The 
disadvantage to studying Shelikof Strait is that it is 
sheltered from the deep ocean swell; this limited the 
wavelength and amplitude of gravity waves available to 
illuminate shoaling features. However, because it is 
representative of likely candidates for deep-water shipping 


Fóltes, 1t is a practical area to demonstrate the 
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feasibility of use of SAR to locate possible hazards to 


navigation in hydrographic presurvey planning. 
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II. SYNTHETIC APERTURE RADAR IMAGERY 


A. SAR IMAGING PRINCIPLES 

A SAR transmits a coherent electromagnetic pulse which 
is reflected back to the satellite by resonant (Bragg) 
scattering at the ocean surface. The Bragg scattering cross 
section, and therefore the strength of the returned pulse, 
depends directly on the spectral energy density (per unit 


area) of ocean surface waves of the Bragg wavelength 
La A 72 sinO 


where A, is the radar's electromagnetic wavelength and © 
is the incidence angle measured from the normal to the sea 
Surface. Stated simply, as the amplitudes of ocean waves of 
the Bragg wavelength grow larger, the strength of the 
reflected radar pulse increases, and vice versa. 

The SAR emitted a 33.4 microsecond pulse at a wavelength 
of 23.5 cm at surface incidence angles spanning 19 to 25 
degrees (Vesecky and Stewart, 1982). Other characteristics 
of this instrument are summarized in Table 2.1. 

From the data in Table 2.1, it is apparent that ocean 
waves of length approximately 30 cm are responsible for 
Bragg scatter return measured by the SEASAT SAR. Ocean 


surface waves of this length are called "ultragravity 


T9 


waves", which refers to a region of the spectrum where 
surface tension and gravity are of comparable importance as 
restoring forces. Variations in the amplitudes of BO wen 
ultragravity waves may result from variations in wind 
stress, interactions with larger waves (a mechanism which 
helps visualize ocean swell in SAR images), surface stress 


variations associated vith 


Table 2.1. Seasat SAR System Characteristics 


(Vesecky and Stewart, 1982) 


Satellite altitude 800km 
Imaged swath 100km 
Angle of incidence 19 = 
Resolution 25-40m 
Radar wavelength ZO 
Contrast ratio 9db 
Power, average 55W 
Chirp pulse length 33.4 s 
Bandwidth 19MHz 
PRE 1463-1640Hz 
Antenna 
Dimensions LOs. AO 
Gain 35dbD 
Beam width IAS = 


convergence and divergence regions in surface flow, current 
wave interactions, and spatial variations in surface tension 
due to surfactants such as organic oil slicks (which are 
ubiquitous on the sea surface). Several of these mechanisms 
are, in turn, affected by variations in bottom topography, 


wave refraction, or both. 
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Cross track resolution is obtained by "range-gating" the 
return signal into successive time bins of length dt, which 
must be larger than the radar's transmitted pulse length. 
Referring to Figure 2.1, the cross-track range resolution in 


a single cell gated over time dt, is given by 


dx = c dt/(2 sinQ) 


where c is the speed of light and the other terms are 
defined in Figure 2.1. For the SEASAT SAR, a cross-track 
resolution of 25 m is associated with a range gate of 
approximately 62.4 ns; a pulse of this effective length was 
achieved by compressing the 33.4 us chirped pulse (Vesecky 
and Stewart, 1982). 

A fine enough resolution to also resolve 25 m in the 
along-track direction is obtained, with a SAR, by 
synthesizing a very long radar aperture through Doppler 
pascessings For an idealized monochromatic (a pure Single 
frequency) radar, the only portion of the returning pulse 
which has the same frequency as the emitted pulse is that 
portion returning from the surface area perpendicular to the 
spacecraft's direction of travel. The portion of the pulse 
returning from the area forward of the 90° azimuth is 
Doppler-shifted to a higher frequency, and the return from 
the area behind 90° azimuth is Doppler-shifted to a lower 


frequency (Fig. 2.2). The Doppler change in frequency is 
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nearly linear with small angles (in radians) measured from 
W. Therefore, discrete frequency and amplitude returns 
for each transmitted pulse can be assigned to a particular 
surface area. The information from a number of returns from 
a Single surface area is then use to synthesize an image as 
if it were collected from a radar having a long aperture. 
ang 2.3). 


The along-track resolution for SEASAT is approximately, 
dy = A, R/(2 L_cos6) » 


where R is the range (approx. 850km) and L is the length of 


the synthetic aperture (Vesecky and Stewart, 1982). For, 
L = Vt = 17km, 


where V is satellite velocity (7.2km/s) and t is the sample 
length (2.3s), the along-track resolution is approximately 
6.6m. However, SEASAT signal data was processed such that 
four resolution cells in the azimuth direction were 
incoherently averaged for a resulting resolution of 25m. The 
reader is referred to Alpers (1983), Tucker (1983) and 
Vesecky and Stewart (1982), for a more thorough discussion 
on the technical details of the SAR imaging process 
(including the extension of the Doppler processing concept, 


uSing phase coherence to sort the frequency spread of the 
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Figure 2.2. Doppler range of radar return signals. 
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Figure 2.3. A Synthetic Aperture Radar and its parts. 
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transmitted radar pulse from that associated with Doppler 


shifts). 


B. OCEANIC PHENOMENA 
It has been empirically determined that L-band 
backscatter from the sea surface iS proportional to 


y9: 520.1 


, Where U is the surface wind velocity (Thompson, et 
al., 1983). If local winds are less than approximately 2 
m/s (Kasischke, 1982), the surface will not have sufficient 
roughness to produce significant Bragg scattering; such 
areas will appear dark in SAR imagery. The reverse is also 
true. If the surface wind speed is greater than 
approximately 10 m/s, the number and magnitude of Bragg 
waves on the surface will saturate, resulting in a strong 
signal return devoid of distinct spatial features within the 
affected area. 

There are a number of oceanic and atmospheric phenomena 
which affect the radar cross-section of Bragg waves. Some 
of these, (such as wind stress, currents, internal waves, 
wave diffraction and frontal boundaries) produce SAR image 
signatures without interacting with bottom topographic 
features (Kasischke, et al., 1982). In this study, we are 
concerned with discriminating phenomena which create 
distinctive SAR image signatures through interaction with 


bathymetry. 
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There are three categories of phenomena by which bottom 
topography may affect the amplitudes of Bragg waves and 
thus, the radar cross-section, seen as brightness 
modulations in a SAR image. These are refraction of ocean 
surface gravity waves propagating across shoaling 
bathymetry, interaction of surface waves with currents 
modulated by flow across bathymetric features (either 
directly seen by interaction of the the current shear with 
the Bragg waves or indirectly through generation of internal 
waves) and bathymetric forcing of flow which produces 
distinctive surface temperature patterns which affect the 
stability of the atmospheric planetary boundary layer. 

l. Gravity Wave Refraction 

A surface gravity wave propagating into shoaler 
water experiences wave refraction as a function of depth. If 
the changes in On and direction can be measured, the 
depth can be calculated (see section III.D). 

Optical Fourier Transforms of SAR images can be used 
to measure the changes in ocean wavelength and direction 
(Shuchman, et al., 1979a, McLeish, et al., 1980, Shuchman 
and Kasischke, 1981 and Meadows, et al., 1983). The wave 
direction derived from the SAR spectra will be shifted from 
the true wave direction as a function of the wave's 
direction relative to the SAR's viewing axis (Pawka, et al., 
1980). However, for this study the change in relative wave 


directions is used to estimate the general orientation of 


ae 


bottom features. Therefore, the direction bias due to wave 
motion should not affect the conclusions. 

Ocean currents can also cause surface gravity wave 
refraction. Hayes (1980) showed the Gulf Stream refracted 
shoreward propagating waves enough so the change could be 
detected by SAR. Therefore, a certain amount of caution 
must be exercised when making estimates of depths from the 
magnitude of wave refraction. 

2. Non-Linear Gravity Wave Interaction Over Shoals 

Sharp changes in the image tone of SAR imagery have 
been related to shoal areas surrounded by deeper water 
(Kasischke, et al., 1980). Kasischke, et al., (1982) 
postulate that these brightness anomalies are associated 
with abrupt shoaling of the waters. It is surmised that the 
abrupt shoaling causes a rapid transition of gravity waves 
from deep water waves to shallow water waves, which in turn 
strains the ultragravity wave field and changes the radar 
backscatter returned to the sensor. 

3. Changes ot neu meme Velocity 

Several SEASAT SAR images appear to contain direct 
images of bottom topography features (Kasischke, et al., 
1983, Valenzuela, et al., 1983). However, electromagnetic 
waves of the length generated by the SEASAT SAR, 23.5cm, 
will only penetrate a few centimeters into the sea surface. 
Therefore, SAR must be imaging the results of an interaction 


between the bottom topography and an oceanic process 
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manifested on the ultragravity wave field on the ocean 
surface. 

Conventional radar assumes either the sensor itself 
SEetene observed object is stationary. If a water sumface 
element of high radar cross-section is held stationary on 
the sea surface and observed by SAR as the satellite crosses 
overhead the Doppler shift in the returning signal phase ls 
uEBunctron only of the azimüth angle: In this case, the 
Doppler shift can be properly taken into account during 
processing (Shuchman, et al., 1981, Vesecky and Stewart, 
1982, and Tomiyasu, 1976). Thus, the return signal from a 
particular water surface element can be determined from its 
Doppler shift and corrected to its instantaneous location in 
the radar beam as shown in Figure 2.4. 

If the entire surface were to move at a constant 
velocity during the time the radar beam was ss dO- 
Doppler shift would be modified shifting the placement of 
the return signal intensity as shown in Figure 2.5 (Vesecky 
and Stewart, 1982). Note that the surface velocity vector 
must have a component in the azimuth direction or the 
Doppler shift will not be altered. 

If the surface element having increased roughness 
were to move at a faster rate than the adjacent elements 
(i.e., a jet), the return signal intensity would again be 
shifted (Fig. 2.6): we would now observe an energy peak 


and adjacent minimum, because the adjacent element is not 
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Figure 2.4. Signal intensity returned from a stationary surface. 
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Figure 2.5. Signal intensity returned from a surface with 
moving elements of constant velocity. 





Figure 2.6. Signal intensie, returned from a surface with moving 
elements of varying velocity. 


filling in the 'space' for that particular Doppler return. 
As seen in SAR images, the energy peaks appear as brighter 
pixels and the minima appear as darker pixels. 

In a SAR image of a region of strong surface current 
shear, velocity bunching tends to produce adjacent bright 
and dark bands across the edge of the shear zone. This 
effect, together with current-wave-bathymetry interactions 
produce SAR images of bathymetric relief. 

Nantucket Shoals, where strong tidal currents flow 
across distinct shoals, has produced some of the more 
dramatic examples of this type of signature (Fig. 2.7). 
Environmental conditions required to observe this phenomenon 
at this site have been determined as: (1) a tidal current 
greater than 0.4 m/s, and (2) a wind velocity between l and 
7.2 M/s (Kasischke, et al., 1983). 

Figure 2.8 Bis a schematic illustration of the 
interactions between tidal currents, bottom topography and 
capillary waves, and is labelled to suggest the resulting 
variation in SAR image intensity. The dark tones of the SAR 
imagery correspond to areas where the current velocity 
increases due to a decrease in water depth, decreasing the 
amplitude of the ultragravity waves. The bright tones occur 
where increasing depth has slowed the current velocity, 
compressing the ultragravity waves (increasing their 
amplitude). The gray tones occur over areas of uniform 


depth, where the wind is able to restore the ultragravity 
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Figure 2.7. SAR imagery of Nantucket Shoals - 
pseudo-imaging of bottom features. 
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Schematic diagram of current, bathymetry and ultra- 


gravity wave interactions and the resultant SAR image intensity 
variation (Kasischke, et al, 1983). 
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waves to equilibrium. The reader is referred to Kasischke, 
et al., (1983), for a more exacting description MUS 
process. 

Kasischke, et al., (1980) determined general 
correlations between this type of signature and the charted 
bathymetry for five areas: (1)DMA-Bahamas calibration area, 
(2) Southern Bahamas, (3) Nantucket Shoals, (4) English 
Channel, and (5) North Rona Rock. The correlation was done 
manually by placing a transparent reproduction of a chart 
over the SAR imagery. If there was a 20 gecuccheme dane meg 
water depth within a few hundred meters of a possible 
feature as shown on the SAR imagery, a correlation was 
assumed to exist; correlations for the five areas ranged 
between 60 and 100 percent. 

Kasischke, et al., (1982) performed a more rigorous 
Study of the relationship between SEASAT SAR T and 
Water depths in Nantucket Shoals. They found a good 
correlation between signal intensity and water depth along 
single range lines, (Fig. 2.9). This result suggests that 
for well-defined situations, it may be possible to 
distinguish relative depth changes between discreet points. 

4. Coastal Upwelling 

In general, water in upwelling areas is colder, less 
saline (U.S. west coast, particularly) and richer in 
nutrients, all of which tend to decrease surface tension 


(Knauss, 1978). The decrease in surface tension, in turn, 
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Figure 2.9. Comparison of SEASAT, LANDSAT and Depth data along a 


range line (Kasischke, et al, 1982). 
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decreases the energy density of capillary and ultragravity 
waves on the surface, and therefore decreases the radar 
cross-section. Hence, ocean upwelling areas will have a dark 
appearance on SAR images. 

Upwelling produces modulations in SAR images, both 
from surface dynamics associated with fronts, and due to 
changes in wind stress associated with the effect of changes 
in air-sea temperature differences (and thus, atmospheric 
planetary boundary layer stability), over cold upwelling 
regimes. Further, coastal upwelling is associated with 
topographic structure of the continental shelf and slope 
(Mooers, et al., 1977), and intense upwelling centers are 
often associated with submerged coastal headlands (Arthur, 
1965) and over the heads of submarine canyons (Shaffer, 
1976). Kasischke, et al. (1983) have related the appearance 
of a frontal boundary on a SAR image (revolution 762, August 
19, 1978) to upwelling produced by current flow over a shelf 
in the Faeroe Islands Thompson and Golding (1981) showed 
that strong tidal suctions around steep bottom slopes may 
induce upwelling. Hence, the appearance of an upwelling 
signature in SAR imagery may suggest the presence of a 
shallow bottom sloping steeply into deep water, or of a 
submerged headland (e.g. Arthur, 1965). 

5. Internal Waves 
The detection of internal waves by SAR has been well 


documented, (Trask and Briscoe, 1983, Alpers and Ettore, 
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19857 Hughes and Gower, 1983, Kasischke, et al;, 1982, 
smucnman and Kasischke, 1979, Beal, et al., 1981 and Allan, 
1983). There are several processes which can generate 
imeernal waves, including wind stress, surface wave 
interaction, baroclinic instability, Ekman layer instability 
and current flow over a bottom feature. However they may be 
generated, internal waves create surface convergence/ 
divergence zones which modulate the Bragg wave field at the 
sea surface, and thus produce patterns of dark and bright 
bands recognizable as internal wave train signatures in SAR 
imagery. 

The shape of a given internal wave signature is 
often dependent upon the shape of the bottom feature causing 
the internal wave. In all cases, the signature is a 
distinct set of wave packets. In the Hudson Canyon Area, 
where the signature Meis over a ridge (Kasischke, et al., 
1980) the wave packets are slightly curved and appear to be 
aligned such that each packet follows one after the other. 
In the present study where an internal wave signature 
appears over a seamount (IV.B.2.), the wave packets are much 
more curved and do not appear lined up in a row. From 
observations made in section IV.B.l., the effects of current 
shear are apparent in the wave packet. 

6. Fresh-water Runoff 
Fresh water runoff from rivers and estuaries is 


often of different temperature than the ambient ocean water. 
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In the Shelikof Strait area we would expect fresh-water 
runoff to be associated with snow melt along the Alaskan 
Penninsula, and we therefore expect to find cold, fresh 
water pools and plumes overlying warm, saline ocean water. 
We therefore also expect a stable atmospheric planetary 
boundary layer, reduced wind stress and low radar 
cross-section over fresh water areas. Conversely, in the 
ambient waters of Shelikof Strait which are warmer and 
saltier, we would expect a less stable planetary boundary 
layer, larger wind stress and therefore, brighter patterns 
in the SAR imagery. 

A light, fresh water plume entering a nearshore area 
over nearly flat topography will tend to propagate offshore, 
but the plume will be arrested at an abrupt shelf break. 
Vorticity conservation across an abrupt deepening of the 
water column will tend to deflect the flow to the right (in 
the northern hemisphere), and establish a front along which 
the flow P proceed parallel to the shelf break in 
approximate geostrophic balance. This phenomenon will tend 
to trap fresh water Te Ore flat shoals, and retain 
distinct areas of fresh water, the boundaries of which 
closely resemble the edge of the shoal. Such phenomena in 
SAR images may be interpreted as indicating a trend of 
steeply deepening bathymetry in the general location of a 


shoal area. 
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III. DATA AND METHODS 


A. TEST DATA 
l. SAR Imagery 

SEASAT acquired the SAR imagery used in this study 
on three orbits over Shelikof Strait, Alaska between July 
17 and August 22, 1978. Details of the system 
characteristics are shown in Table 2.1. The optically 
processed imagery of interest was obtained from the 
Environmental Data Information Service, National Climate 
Center, Room 100, World Weather Building, Washington, DC 
PASS, tel. (301)763-8111. 

Shelikof Strait was chosen as the study area due to 
the availability of SAR images of the site from three 
separate SEASAT revolutions and because it is typical of 
possibly useful deep-water shipping routes. The imagery 
from revolution 289 on July 17 was captured on a descending 
path with an azimuth of 212 degrees. The site was imaged 
again on ascending passes on July 19 (revolution 323) with 
an azimuth of 328 degrees, and on August 22 (revolution 811) 
with an azimuth of 328 degrees (Fig. 1.1). 

p xtompEoln Data 
A secondary factor in selecting Shelikof Strait as a 


study area was the availability of bathymetric maps of the 
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area.  Bathymetric maps are a better data source than 
navigation charts for this type of study. A standard 
navigation chart depicts discrete depths which are 
representative of the soundings taken in that area. The 
depth represented on the chart is always the shoalest in the 
area, So that the actual depth for a given point may be 
somewhat greater than that shown. Bathymetric maps, on the 
other hand, are contoured directly from survey soundings (as 
are contour lines on topographic maps) and present a more 
accurate description of the bottom topography. 

The control data for this study were bathymetric 
maps and hydrographic charts obtained from the National 
Ocean Survey, National Oceanic and Atmospheric 
Administration, Rockville, MD 20852. The three bathymetric 
maps covering the study area inside Shelikof Strait are 
Afognak (NOS No. 5-4), Iliamna (NOS No. 5-1) and Mt. Katmai 
(NOS No. 5-3), all at 1:250, 000 "scale: Historica E e 
these charts show hydrographic surveys covering selected 
subareas were performed within the past 20 years at survey 
scale 1:40,000 or greater, with line spacings generally less 
than 0.5 nautical miles (Appendix B). 

The hydrographic chart used for study area VIII, 
southeast of Kodiak Island, is a small-scale (122,100,000) 
hydrographic chart NOS 531, titled 'Strait of Juan de Fuca 
to Kodiak Island,' tenth edition, printed in 1975. The 


hydrographic charts used for Nantucket Shoals were: NOS NO. 
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13237, 'Nantucket Sound and Approaches,' scale 1:80,000, 
28th edition printed in 1981, and; NOS No. 13241, 'Nantucket 
iewand, scale 1:40,000, llth edition published 1981, 

3. Oceanographic and Meteorological Data 

Oceanographic and meteorological phenomena that are 
relevant to this study include ocean currents, tides and 
tidal currents, winds, and precipitation . Some knowledge 
pEEneSse conditions prevailing during each observation, sand 
an understanding of how they can affect the Bragg wave radar 
cross-section, is essential for ascribing SAR signatures to 
the proper oceanic process. 

The mean current flow through Shelikof Strait was 
from northeast to southwest during late July and August 1978 
when SEASAT acquired SAR imagery of the area. The 
along-shore mean current speeds ranged from 10.2 to 17.6cm/s 
in the upper 25 meters of the water column.  (Mysak, et al., 
198l and Muench, et al., 1978) 

Muench and Schumacher (1980) constructed an 
empirical cotidal chart for the principal lunar diurnal (M5) 
Ed constituent for Kodiak Island and Shelikof Strait area 
(Fig. 3.1). The chart was constructed by interpolation 
where tide stations were closely space relative to the tidal 
wavelength and by estimation where tidal observations were 
scarce. This chart is used here in conjunction with 
published tidal charts for Kodiak Island to estimate the 


tide stage associated with each SAR image. 
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Figure 3.1.  Cotidalichart for M tide'component. Solid lines 
are cophase lines referenced to Green A Dashed lines are 
coamplitude in centimeters (Muench and Schumacher, 1980). 
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Tree are no Specific windmor preeipitation data 
avatlable for Shelikof Strait during the tim@@ of SEASAT 
coverage. Some information can, however, be inferred from 
inspection of the images. 

Localized rain squalls are more likely to affect 
bathymetric interpretations than a generally widespread rain 
storm by creating local contrasts in the radar backscatter. 
However, rain squalls have a very distinct SAR signature 
(Kasischke, et al., 1982) and the collected images do not 
have any of these signatures. Therefore, precipitation 
Should not affect the bathymetric interpretations of this 
Study. 

Land topographic features affect wind stress 
patterns, which may in turn affect the spatial patterns of 
the Bragg wave radar cross-section. There are a number of 
along-shore areas in the Strait and a number of bays and 
inlets which had a low backscatter intensity in the imagery. 
Apparently these areas were protected from the wind by 
topographic features. These areas were excluded from the 
analysis because of the low radar return and lack of signals 


Weerul for inferring bathymetric information. 


B. STUDY AREAS 
The eight subareas used for analysis in this effort are 
shown in Figure l.l. The subareas were selected for 


analysis on the basis of two criteria: 1) observed visual 


45 


signatures in the SAR imagery, and/or 2) the likelihood of 
finding significant wave refraction (as estimated from 
bathymetric maps). 

The observed visual signatures in the imagery were 
transferred to a MYLAR overlay. The overlay was then 
superimposed on a bathymetric map using a zoom transfer 
scope by colocating topographic features (e.g. shorelines, 
islands, inlets). In this way, geometric distortions in the 
SAR imagery were corrected for and the signatures could be 
correlated to bathymetric features. 

The wave refraction data was collected in a similar way. 
Once general areas of expected refraction were identified, 
OFT's of the areas were measured and the area of each OFT 
was transferred to a MYLAR overlay. M Each OFT area was then 
colocated with the proper bathymetry using the zoom transfer 
scope and bathymetric maps. 

In general, the areas chosen for OFT analysis had a 
nearly linear bottom slope. in Some casos SUCIA SC 
the slope was fairly shallow. Area III, on the other hand, 
had a fairly steep slope, and the other areas fell between 
these two extremes. 

There were two exceptions to this criteria. Area II 
includes a large seamount in the northeast entrance to the 
Strait; an attempt was made to measure the surface wave 
structure as it passed over the seamount, with some success 


(IV.B.2). In Area VIII, a large-scale natural channel 
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approaches Kodiak Island from the southeast. An attempt was 
made to locate this channel by the change in surface 
wavelengths as measured using the OFT, again with modest 
success (IV.B.8). 

Originally, there were 15 areas chosen for OFT analysis. 
However, in seven of these areas the OFT did not yield 
discernible information on the surface wave structure, and 


therefore they were excluded from the analysis. 


C. VISUAL INSPECTION OF FEATURES 

l. Internal Wave aces 

Any modulation or disturbance of the Bragg wave 

amplitude on the sea surface will be seen in the SAR imagery 
as a brightness modulation or anomaly. Previous studies 
have shown that major subsurface features, such as seamounts 
and canyons on shelf edges, often are associated with 
generation of internal waves with surface signatures 
consisting of dark and bright bands (Shuchman and Kasischke, 
1979). These signatures are manifested in the imagery as 
wave packets radiating outward from the structure that 
caused the deflection. Internal waves which systematically 
modulate the surface radar cross-section are often generated 
Dn ternal tildes (Apel, et al., 1975) and other barotropic 
currents (Kasischke, et al., 1983) interacting with shelf 


breaks and other bump-like features in the bathymetry. 


47 


ZAS 

A surface wave structure incident upon an opposing 
surface current will be refracted against the current 
velocity vector, and if the current is flowing in the 
direction of wave propagation the waves will be pushed 
towards the normal to the current velocity vector (Fig. 
3.2). Hence, the appearance of an apparently sheared wave 
structure can indicate the presence of a current. 
Furthermore, internal wave packet signatures extending 
across regions of current Shear will also be distorted 
through differential advection. 

3. Non-Linear Gravity Wave Interaction 

Gravity waves propagating across broken, rough shoal 
areas experience non-linear interaction with the bottom. 
This non-linear interaction could be caused by waves 
crossing rocks which are classified on navigation charts as 
‘rocks awash' or ‘uncovered’ (at low tide). These 
interactions would cause the surface gravity waves both to 
break and to undergo a rapid transition from a deep-water 
wave to a shallow water wave. Such a rapid transition tends 
to strain the ultragravity wave structure, and consequently 
modulates the radar cross-section, producing a much higher 
radar backscatter over the shoal than in the surrounding 


area (Kasischke, et al., 1983). 
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Figure 3.2. Refraction of wave structures incident upon a 
current (Phillips, 1981). 
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D. OFT ANALYSIS OF SURFACE GRAVITY WAVES 
l. Wave Refraction 

Linear Wave Theory is pm In this study +t o predi È 
water depths and shoaling using measurements obtained from 
OFT of SAR imagery. Linear Wave Theory in shallow water 
predicts that as waves propagate shoreward, wavelength 
decreases with decreasing water depth. 

For shallow water waves where the ratio of water 
depth to wavelength is less than 1/25, the water depth is 


related to wavelength by the equation 


where d is water depth, L is wavelength, T is wave period 
and g is the gravitational constant (Coastal Engineering 
Research Center, 1984). For deep water waves, where the 
ratio of depth to wavelength is greater than l/2, the wave 


period T is given as, 

SIE T = (2gL.)/g, 

where Lo is the deep water wavelength (Coastal Engineering 
Research Center, 1984). The wave period remains constant, 


so we may substitute equation 3.2 into equation 3.1 to 


obtain, 
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Since L and L can be measured from OFT of SAR images, it 
may be possible to determine local water depths using 
equation 3.3. 

The change in the direction of wave propagation due 
to refraction can also be used to determine water depth. As 
waves propagate shoreward at some angle to shoaling bottom 
contours, they tend to turn in the direction of the 
shoreline as shown in Figure 3.3. Snell's Law relates the 


direction of propagation to the wave velocity as 
3.4 i n o< inc = 
sin </sin = C M Cer 


where oO, c, are deep-water direction and velocity, 
respectively and Er C, are those at isobath depth Za The 
change in phase speeds is related to water depth and 


wavelength as, 
$55 cmon an ( 2 Ls), 


(Coastal Engineering Research Center, 1984). 
Memmay thus combine 3.4 and 3.5 to relate water depth amd 


wavelength to propagation direction as 


3526 sino/sinos = tanh( (24d)/L). 


al 





Figure 3.3. Refraction of wave structures incident upon 
bathymetry. 
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The left-hand ratio of equation 3.6 and wavelength 
can be obtained from OFT's of SAR images. Using that ratio 
and observed wavelength, water depth can be determined from 
empirical tables published in Shore Protection Manual 
(Coastal Engineering Research Center, 1984). 

Equation 3.1 may also be used to measure relative 


shoaling. For any discrete shallow water area, 
E 2 2 

3m 7 di = Li AT is 

And for an adjacent shoreward area, 
T 2 2 

Sg dj 7 L, INTA 


Dividing equation 3.8 by equation 3.7 yields, 


A 2 2 
ory d/d} = Lo /Li ; 


where Li and L^ are different wavelengths of the same wave 


group for adjacent areas. 


E Optical Fourier IBAS cons 


Coherent, collimated light passing first through an 
object transparency and then through a transform lens 
separates the object pattern into its Fourier components at 


the transform lens's focal plane (Goodman, 1968). The 
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Fourier components of a SAR image are thus displayed in 
wavenumber (k) space and the intensity of the transform 
image at wavenumber (kj ko) is proportional to the squared 
gray level amplitude of features at that vector wavenumber 
in the image. 

The distance from the transform origin (at the focal 
plane) to a component 1s thek Spacewaye mumk NE 
particular component. Also, the directional orientation of 
a component around the origin is determined by the wave's 
direction relative to the SAR sensor. Shuchman and 
Kasischke (1981) used these principles to predict water 
depth from OFT spectra of SAR imagery. A detailed 
description of the optical transform procedure is given here 
in Appendix A. 

A helium-neon laser was used to perform the OFT's by 
focusing the beam through one lens onto a 1 em area of an 
image transparency, then through a second lens. The OFT 
appeared on the focal plane of the second lens and was 
recorded on 35 mm infrared film. 

The K-space wavenumber was chenge asure OM 
mathematically transformed to the corresponding wavelength 
for the surface gravity wave. The angle of the component 
was also measured (with respect to the satellite's heading) 
and converted to a wave direction. 

DHecerror assocratedwT tH LDEMRLCUIEOCCLE US 


was on the order of +l mm. A geometric correction was made 
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Ae canstorm distortion due to the camera offset. A 
Sempesite®eorrection for errors! aseec iat cae@with tite Sept ileal 
train was determined empirically (Appendix A) and applied to 


eaem OFT results 


E. ASSESSMENT OF DATA QUALITY AND ERROR SOURCES 

The imaging medium itself contains inherent errors. The 
original data was optically processed and printed on 
photographic film, and then copies were generated and sent 
to the Environmental Data Information Service (EDIS) for 
archival. When EDIS receives a request for imagery, they 
generate copies from the archived film and the user receives 
a third-generation film product. Moreover, in the original 
printing of the imagery the feed mechanism which ran the 
film strip through the print operation varied in its speed, 
creating a randomly varying distortion in the azimuth 
Gateection. 

Another error source lies in the dynamic capability of 
the Heu The dynamic range of the SEASAT SAR sensor 
was 60db. However, the photographic film's range is 17-20db. 
Therefore, the original data has been compressed with an 
ensuing loss of information (Kasischke, et al., 1980). 

There are two other inconsistencies in the film medium. 
he first is an artificial darkening of the gray tone from 
one side of a subswath to the other (in the range 


direction). The second is a change in the scale of the 
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imagery, from 1:480,000 at near-range to 1:520,000 at 
far-range. However, both inconsistencies are insignificant 
at the local scale and have negligible effect on the results 
of this study as only a small portion of each subswath was 
used for each evaluation of a bathymetrically related 
signal. 

The use of digitally-processed SAR images would not have 
significantly improved the accuracy of the quantitative and 
signature analysis using the optical techniques employed in 
this study. Even if the images were digitally processed, 
the same dynamic range restrictions would apply when the 
data were recorded on film to create the image. Geometric 
distortions in the azimuth direction would not be present; 
however, they were adequately compensated for in the 
optically-processed imagery. The other two T 
noted, artificial change in tone and a changing scale, would 
not have been present, but as already noted, these effects 
were not a significant error source in the present study. On 
the other hand, alternative approaches using interactive | 
digital image processing systems could take advantage of the 
improved dynamic range and signal-to-noise ratio inherent in 
the digital data. 

The survey line spacings from which the bathymetric maps 
were compiled were generally 0.1-1.0 nautical miles 
(Appendix B), which transforms into 0.4-4.0 mm spacings at 


image scale. The area used for the OFT process was 1 cm; the 
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wave structure detected, therefore, could not be located 
more precisely than within that one square centiméter. 
Therefore, the survey line spacings should not have affected 
the OFT results. 

Tne application Of Optical Fourier Tramsform theory 
makes several assumptions: 1) that the lenses are perfectly 
ground, 2) that the laser beam is perfectly collimated, 
producing perfectly parallel rays, 3) that the focusing 
aperture is perfectly square, and 4) that the transparent 
portion of the transparency has no thickness. In this 
study, all of these assumptions were incorrect. In order to 
empirically compensate for these errors, OFT'S were 
performed on gratings with known spacings and the predicted 
wavenumbers were compared to the measured wavenumber. A 
correction coefficient was then calculated and applied to 
all measured wavenumbers when calculating the corresponding 
wavelength. This coefficient did not correct for any 
distortion of the wave direction, however, which may have 
contributed to the relatively poor performance of wave 
direction changes as a depth predictor. 

There were possible colocation errors associated with 
transferring points from the SAR image to the corresponding 
bathymetric map. The SAR image is a slant plane 
representation of a curved, three-dimensional surface, with 
a changing scale from near range to far range. The 


bathymetric map is plan view representation of a curved, 
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three-dimensional surface with a changing latitudinal scale. 
These errors were minimized by using a zoom transfer scope 
to perform an analog "rubber sheet" stretching to colocate 
identifiable features such as shorelines, swampland and. 
islands. Because the water areas studied here are 

'" surrounded by well-defined landmarks (except for Area VIII), 
the colocation error was estimated to be «5$ of a given OFT. 
The uncertainty may be as high as 10% of an OFT area in Area 
VIII, which is outside the Strait. In no case is this a 


Significant source of error. 
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IV. ANALYSIS 


The SAR images, bathymetric data and oeanographic 
information described in Chapter 3 are analyzed in this 
chapter for subareas within Shelikof Strait. The case study 
subareas were EINE according to the criteria outlined in 
Section III.A. This chapter will first discuss the physical 
environment (including bathymetry, oceanography and wind 
conditions) in Shelikof Strait at the times the SEASAT SAR 
images were acquired, and then data for the eight subareas 


are analyzed on a case-by-case basis. 


A. OVERVIEW-OF THE PHYSICAL ENVIRONMENT 


. 


1. Geographic Setting and Bathymetry 


In this section, the overall bathymetric topology of 
the Shelikof Strait area is described, followed by a brief 
discussion of the relatively sparse information available on 
tides, tidal currents, mean current, and density 
stratification of the water column (including the possible 
influence of fresh water runoff). Finally, the prevailing 
wind conditions during each of the image acquisition periods 
are discussed, based both on synoptic scale surface wind and 
pressure analysis, and on visual inspection of the overall 


SAR images. 
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Shelikof Strait is located between Kodiak and 
Afognak Islands, and the southern coast of the Alaskan 
Peninsula, with a northeast-southwest orientation. The 
principal bathymetric features are shown in Figure 4.1. 

The Strait has a depth of 196 m in the northeast 
entrance, increasing to 340 m in the southwest. The Strait 
is about 33 km wide with a generally smooth bottom, a gentle 
downward slope to the southeast in the upper portion of the 
Strait, and a slightly steeper upward slope to the southeast 
in the lower portion. The sides of the Strait are generally 
steep, shoaling from 180 to 50 m over a distance of less 
than 1.8 km, creating a coastal shelf with steep sides. The 
top of this shelf has a gentle slope and extends 5 km out 
from the coastline on the southeastern side of the Strait 
A UD km on the nort bel Side. The mao cs mane 
feature in the strait is a seamount which shoals from 160 to 
100 m located in the northeastern entrance to the Strait. 

2. Tides and Currents 

Tidal Eon used in this study is derived from 
tide tables published by the National ocean Survey for 19787 
and from the cotidal chart in Figure 3.1 (Muench and 
Schumacher, 1980). All extrapolations are referred to the 
predicted high and low tides for the tidal station at the 
town of Kodiak, Alaska (Figure 1.1). The tidal stage for 


each image acquisition period is presented here, but no 
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Figure 4. 
Alaska. 
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l. Principal bathymetric features in Shelikof Strait, 
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attempt is made to generate tidal heights, due to lack of 
reliable tidal data inside Shelikof Strait. 

The mean current in Shelikof Strait flows to the 
southwest, and is driven by two physical phenomena: a 
pressure gradient set up by the strong southwesterly flowing 
Alaskan Stream as it passes to the south of Kodiak Island, 
and a southwesterly continuation of the Kenai Current, which 
originates east of Kennedy and Stevenson Entrances (Fig. 
4.1) (Mysak, et al., 1981). Seasonal variations in the 
Kenai Current cause, in turn, seasonal bd omncm Me 
mean current flow through the Strait. The mean velocity 
through the Strait at the 25-meter depth was 10.2 cm/s in 
early summer 1978 and 17.6 cm/s in late summer 1978 (Mysak, 
et al., 1981). Further, there was a 6 cm/s current velocity 
associated with the M, tidal component superimposed on the 
mean flow. (Mysak, et al., 1981) 

Normally, at low tide the tidal current has a 
reversed direction with similar velocity compared to tidal 
currents during high tide, and slack water occurs 
approximately halfway between high tide and low tide 
(Redfield, 1980). However, when the tidal current is 
superimposed on a consistent current flow, such as in 
Shelikof Strait, it is difficult to estimate the effects of 
tidal currents on net water flow during low tide, 
particularly in areas where the flow may interact with 


bathymetry. This effect may induce a 50% fluctuation in the 
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mean current condition during the season studied here 
(Mysak, et al., 1981). 

The SEASAT SAR imaged the area shown in Figure 4.8 
Hung revolution 289 at 0545 GMT on July 17,1978.  Thd tide 
emes predicted low tides for Kodiak at 0428 GMT Ána 1602 
GMT, and high tides at 1051 GMT and 2226 GMT. It is assumed 
that the differences between successive high and low waters 
Hewequivalent to a 90° change in phase, and that the phase 
difference between Kodiak and Shelikof Strait is 20° (Figure 
3.1). It is estimated that the waters inside the Strait 
were just approaching low tide as Seasat passed overhead, 
and that the mean current flow was ~10 cm/s to the 
southwest. 

SEASAT acquired the SAR images shown in Figures 4.4 
and 4.16 at 1430 GMT on July 19, 1978 during revolution 323. 
The tide predictions for Kodiak were low tides at 0610 GMT 
and 1757 GMT and a single high tide at 1233 GMT. Thus, 
Shelikof Strait is estimated to have been at high tide with 
a mean current flow of ~16 cm/s to the southwest at the time 
the image was zed on revolution 323; 

The third image was acquired at 1700 GMT on August 
SS during SEASAT revolution 8ll (Figs. 4.2, 4.4, 4.7, 
4.14 and 4.19). Predicted tides for Kodiak were low tides 
at 0919 GMT and 2156 GMT and high tides at 0315 GMT and 1539 
GMT, indicating the tide inside the strait was just entering 


its ebb stage during this acquisition period. Again, based 
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on the results of Mysak, et al (T98) thes cubic ow 
estimated to be ^18 cm/s to the southwest. 
3. Density Stratification 

Thermal infrared satellite images of Shelikof Strait 
distinguished a band of surface water that was colder than 
the apparent ambient waters in the rest of the strait 
(Mysak, et al., 1981). This band of colder water appeared 
along the northwestern Shoreline of the strait and appeared 
wider in October than in January. Thus, the appearance of 
colder surface layers may be due to the influence of fresh 
water runoff, which is cold relative to the warmer, saltier 
ambient waters in the strait. 

Although hydrological data for the Alaskan Peninsula 
are not available, we expect snow melt in this area to have 
produced a significant amount of runoff during the summer 
months of this study. This cold, fresh runoff would enter 
the Shelikof Strait from inlets along the coast. Fresh 
water runoff can create stationary ocean fronts over shallow 
water that is bounded by a sharp increase in water depth. A 
cold, cm surface layer inshore of the front would produce 
a local increase in the atmospheric planetary boundary layer 
Stability, so that wind stress would be greater over the 
ambient surface water than over the cold fresh water area. 
This would cause an abrupt change in the radar cross-section 
and hence, the ambient waters would appear brighter in SAR 


imagery than would the colder, less saline waters. Some 
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examples consistent with this paradigm occur in the images 
studied here. 
4. Wind Conditions 

There is no localfyind deta tavanMable for the areas 
of interest within Shelikof Strait. Therefore, wind 
conditions are estimated from National Meteorological Center 
(NMC) Final Analysis data for the SEASAT acquisition periods 
used in this study. 

The data for July 17, 1978 (SEASAT revolution 289) 
indicate winds were moderate from the northeast. Thus, the 
wind stress was roughly aligned along the axis of Shelikof 
Strait. Therefore, topographic effects on wind patterns 
over the strait should have been minimal. 

The data for July 19 and August 22, 1978 (SEASAT 
revolutions 323 and 811, respectively) indicate winds were 
moderate from the southeast, a ENERO mn. The wind 
patterns were most certainly altered by topography as they 
crossed Kodiak and Afognak Islands and as they approached 


the mountainous topography of the Alaskan Peninsula. 


BeeeecAoh STUDIES OF LOCAL SUBAREAS 

The SAR imagery contains organized patterns of radar 
cross-section returns throughout the study area. Four 
subareas in and around Shelikof Strait were chosen for 
analysis based on macro-scale radar brightness features in 


the SAR imagery. Internal wave signatures are visible in 
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the images of subareas I, II and VII. Also, brightness 
anomalies are also visible in subareas II, III, IV, VI and 
VT TS 

Four other subareas, Areas III, IV, V and VIII, were. 
selected for OFT analysis based on the expected wave 
refraction as determined from bathymetric maps of the area. 
OFT analysis was also done on the subareas chosen from 
visual inspection of the imagery. The location of the eight 
case study subareas are indicated in Figures l.l and 4.1. 

l. Area I: Internal Waves and wouimiacer Gna eee 

Refraction 

Area I lies between Shelikof Strait and the mouth of 
Cook Inlet, (Figure 1.1). The imagery analyzed for this 
Subarea was acquired during revolution 811 just before high 
tide, (1700 GMT July 19, 1978). There were moderate, 
southeasterly winds in the area TRU acquisition. 

The transform areas are not close to shore and the 
internal wave signature approaches Augustine Island from 
approximately the same direction as the wind. Coastal 
topography should not have caused large variations in the 
down-strait wind pattern, and therefore, should not have 
strongly influenced the patterns of wind-induced radar 
cross-section in the SAR imagery. 

Figure 4.2 shows an internal wave signature in the 
form of a wave packet propagating toward Augustine Island, 


beginning approximately ten km offshore with a heading of 
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Figure 4.2. Imagery from revolution 811 covering Area I, 
showing transform areas Bl4, Bl15 and Bl6. 
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3209 with wavelengths decreasing from 1.7 to 1.0 km. There 
is a slight northwest shear in the wave structure, which is 
an indication of an incoming current. As the tide is still 
approaching high water, the internal waves are possibly 
incident upon the incoming eka BA LAOS) foe), 

The internal wave packet appears in water with a 
depth of 45 m where the bottom slope is very small (almost a 
plateau) and continues to the island's shelf. The wave 
packet is indicated by curved parallel lines in Figure 4.3. 
Going seaward, the bottom's downward slope begins to steepen 
at around the 60 m curve and increases drastically at the 75 
m curve (Fig. 4.3). The steeper slope may be reflective and 
turning the internal waves back toward deeper water and 
down, with no interaction with the surface and no SAR 
signature.  Shoreward of the 60 m curve the slope is very 
gradual, and may be horizontally transmissive. 

The internal wave signature wavelength decreases and 
brightness contrast increases as the wave approaches the 
more steeply sloping bottom near Augustine Island (Fig. 
4.3). These observed characteristics of the wavetrain, 
e.g., change in wavelength and increased amplitude, are 
consistent with theory (Wunsch, 1968). Also, the observed 
curvature in refraction is consistent with curvature in the 
bathymetry. 

The bathymetric maps show a large number of rocks 


awash near and to the south of Augustine Island, as well as 
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Figure 4.3. Bathymetry of Area 1, (depths in meters, scale 1:250,000). 


a rough, broken bottom to the east of the island. These 
features should have produced distinct signatures in the 
Form or brightness anomalies in the imagery. However, only 
one signature appeared; a small, elongated patch east of the 
island. Tidal heights in this area are typically tens of 
meters in magnitude and the water level was approaching high 
tide during the acquisition period. Hence, lack of 
signatures in the SAR image may be due to too much water 
over these features, preventing detectable modulation of the 
surface radar cross-section. This possibility clearly 
suggests the advantages of obtaining several SAR images of a 
study area at different tidal stages in future research and 
applications. Whatever the reason, there was no indication 
of these known shoals in the SAR image. This negative 
finding highlights the one-sided nature of SAR information; 
it can be trusted to indicate the possible (even probable) 
presence of shoals, but it can not be trusted to 
definitively indicate the absence of shoals. 

Ten OFT's were generated for PEE in this 
subarea. However, only three yielded measurable information, 
indicating that if there was a dominant wave structure in 
most transform areas, it had no components discernible to 
the OFT process. 

The three transforms which showed a discernible 
dominant wavenumber (Bl14, Bl15, and Bl16), are in the areas 


delineated in Figure 4.3, and the measured wavelengths and 
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directions are given in Table 4.1. The spatial relation of 
the transforms and the large differences in measured 
wavelengths indicate the OFT's measured a different wave 
structure in each transform area. Even though a dominant 
wavelength was discernible in each OFT of Table 4.1, all 
evidence suggests that gravity waves in this area were of 
extremely low amplitude and are therefore unsuited for 


detection of wave refraction effects. 


Table 4.1% Wavelengths and Directions in Area I 


Transform Area Lm) ec (degrees) 
Bl4 a 89 348+0.5 
B15 11:925 052419 18+2.0 
B16 9385 TOS 208150 


D Area IN Internal Waves, pnuBness Anomalies and 
Surface Gravity Wave Refraction 


Area II is located in the northeastern entrance of 
Shelikof Strait. Figure 4.4 shows the imagery acquired 
during SEASAT revolution 811 covering this area and Figures 
4.5 and 4.6 show the bathymetry. 

Synoptic wind analysis shows the winds during this 
acquisition period (August 22) were moderate and from the 
southeast, coming into the subarea across Afognak Island and 
through Kennedy and Stevenson Entrances. te does not 
appeal to be any low radar return areas along “the shoreline 


at the southeastern edge of Figure 4.4, which would be 
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Figure 4.4. Imagery from revolut fond li c OVCcPIMOM Ecc cO 
III and IV, showing transform areas A4, A5, A8, A22, A23, 
A26, A27, A30 and A31. 
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Bathymetry of Area II (depths in meters, scale 1:250,000). 


Figure 4.5. 
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Bathymetry of Area III, (depths in meters, 


scale 1:250,000). 


Figure 4.6. 
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indicative of a lee area, suggesting that local wind 
patterns do not make an Important contribution to the 
features seen in this image. 

There are internal wave signatures in the form of 
wave packets in the imagery which correspond to a seamount, 
seen in the bathymetry near 58%45'N, 153°W. The depth in 
the area of the signatures is approximately 210 m, shoaling 
to 70 m over the seamount. The spatial offset of the 
signatures from the seamount resembles offsets found in 
similar circumstances by Shuchman and Kasischke, (1979). 

Mysak, et al., (1981) have shown that the density 
distribution within Shelikof Strait can be approximated with 
a two-layer model, where the upper layer varies in depth 
from 25 to 90 m. The internal wavelength as measured from 
its SAR signature in Area II was 900 m. The height for the 
upper layer (hy) is 90 m and for the bottom layer (ho), 80 
m. The wavelength (L) is much larger than the heights, 
hence, the wave group velocity, 


(€, - Cotgno 
C,coth(kh,) + Y, coth(kh;) 
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can be reduced to , 
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where 
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Using equation 4.1 and the water layer heights determined 


4.1 C = 0.002 


= 


from Figure 4.5, the internal wave group velocity is, 
C = 25 77km/ iae 

The leading wave is 12 km from the top of the 
seamount shown in Figure 4.5. Assuming the internal wave 
generated at the seamount through interaction with the 
incoming tide, the internal wave required five hours to 
reach its present position; this is consistent with the 
estimated tidal phase. 

In addition to the internal wave signatures, the SAR 
image displays a brightness anomaly directly over the 
seamount. Because the seamount is in the center of the 
strait, it is unlikely that this pattern could result from 
wind variance. Therefore, the increase in radar 
cross-section must be associated with some current-wave 
interaction over the shoal. One possible explanation is an 
interaction between the strong prevailing southwest current 
and. the seamount, which is not a single conical shape, but 
two cone-like structures rising 70 m above and to either 
Side of a small canyon. Also, the slope is much steeper on 


the down-current side of the seamount than on the up-current 
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side. This bathymetry (and its interaction with the 
current) may have modified the surface radar cross-section, 
increasing the backscatter intensity seen in the SAR image. 
The brightness anomaly seen in figure 4.5 correlates 

nicely with the location of the seamount. Unfortunately, 
the imagery acquired during SEASAT revolution 289, which 
covered the same area, had such high levels of backscatter 
intensity throughout the area that no information can be 


Gumeained from sit (i.e., the area was saturated). 


Table 4.2.  Wavelengths and Directions in Area II 


Transform Area L(m) OX (degrees) 
A22 a. 180+0.0 
A23 1,160+89 LS RU. 
1,024+76 SUSO 
A26 9564+51 1701,0 


Three OFT were generated for the surface wave 
structure over the seamount. The locations of the 
transforms (A22, A23 and A26), are shown in Figure 4.5 and 
Table 4.2 shows the measured wavelengths and directions. 

Transform A23 has two measurable components which 
may indicate a wave structure undergoing refraction as it 
propagates over a shoaling bottom; the wavelength shortening 
and the direction change is consistent with the bathymetry 
shown. However, the error associated with measuring 
wavelengths of this magnitude makes it difficult to compare 


similar wavelengths, the one mm standard error of 


p 


measurement affects small wavenumbers much more than large 
wavenumbers. Hence, the best that can be predicted for 
transform area A23 from OFT analysis is shoaler water in the 
westernmost portion of the area. 

There are two cases where depth range and a shoaling 
ratio can be predicted; OFT data from A23 and A22 were used 
to predict depths and shoaling ratios in A22, and data from 
A26 and A22 were also used to predict depths and shoaling 
ratios in A22. The ratio of d/L here is larger than 1/25. 
Therefore, the area should be classified as transitional 
water and the following equation from the Coastal 


Engineering Research Center (1984) used, 
Ae L = (gT^)/(2q) (tanh(2gyd/L)). 


From the bathymetry, the depths in area A23 range 
from 125 to 160 m with an effective median depth of 145 m. 
Using equation 4.2 and the measured wavelength from A23, 
L5471,160 m, the wave period (T) is calculated to be 36 sec 
for a depth of 125 m, and a depth of 62 m is predicted for 
area A22. Using the same data, but a depth of 160 m in A23, 
the wave period (T) is calculated to be 33 sec, predicting a 
depth of 76 m in area A22, giving a predicted depth range of 
62 to 76 m in area A22. A predicted effective median depth 
of 70 m was also calculated using the effective median depth 


and the measured wavelengths for A23 in equation 4.2. The 
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actual depth range in area A22 as shown in the bathymetric 
map is 110 to 160 m with an effective median depth of 135 m, 
a 48-percent error in the predicted median depth. 
Similarly, using OFT data from A26, a predicted depth range 
of 94 to 148 m and an effective median depth of 144 m was 
calculated for A22, a 7-percent error in the predicted 
median depth. The apparent reason for the improvement in 
accuracy is that the path from A26 to A22 follows a ray 
path, A23 to A22 parallels the apparent ray path. 

3. Area III: Brightness Anomalies and Surface Gravity 

Wave Refraction 

Area III is located south of the seamount discussed 
above, and Figure 4.4 shows the imagery of this area 
acquired during Revolution 8ll. Figure 4.6 shows the 
bathymetry for this area and the wavelengths and directions 
for the three transforms developed from the imagery (A4, A5 
and A8) are listed in Table 4.3. Note that all transforms 
have more than one distinct wavenumber component, but Figure 
4.6 shows only one arrow for clarity. 

There are some scalloped bands of brightness off 
Raspberry Island (labelled A in Figure 4.6) and off Afognak 
Island (labelled B in Figure 4.6) in the lower portion of 
the imagery shown in Figure 4.4. These anomalies may be 
interpreted as indicators of the escarpment shown in the 
bathymetric map, even though feature A is 2.5 km and feature 


B is 10.0 km from the escarpment. Although the mechanisms 
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which produce the anomalies are unclear, that uncertainty 
would not preclude taking advantage of the fact that the 
shapes of anomalies are indicative of the trend of nearby 


bathymetry. 


Table 4.3.  Wavelengths and Directions in Area III 


Transform, Area Lm) &< (degrees) 

A4 ENDING 184+0.0 
l,116+91 186+0.5 

PSOE 181+0.0 

AS O 174+0.5 
1,342+124 175+1.0 

857+49 Os e 

A8 27289 25 190+1.0 
2,362+413 189+1.0 


Wave refraction data obtained from OFT's were used 
to predict a water depth range for transform area A4 using 
the previously described method. The depth ranges in A5 and 
A4 are 170 to 200 m and 50 to 200 m, respectively and the 
wavelengths were 1,342 m and 1,150 m, respectively. The 
effective median depths for A5 and A4 are 190 m and 150 m, 
respectively. Thus, the predicted depth range in area A4 is 
117 to 135 m with a predicted effective median depth of 130 
m, an error of l3-percent. 

4, Area IV: Brightness Anomallecskan e Urra cN lu, 

Wave Refraction 

Area IV is in the northern corner of the upper 


entrance to Shelikof Strait, (Fig. 1.1). The bathymetry and 
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Erenstorm areas are shown in Figure 4.7. Figure 4.4 shows a 
large, bright area in the imagery to the southwest of the 
three transform areas, which corresponds to an extremely 
shoal area with rocks awash, (Fig. 4.7). The imagery shows 
a fairly rough surface, possibly caused by water crossing 
the submerged rocks at high tide. 

Table 4.4 lists the wavelength and directions 
measured from the OFT's in transform areas A27, A30 and A3l. 
Note that area A3l has two wave components, but Figure 4.7 


shows only one arrow for graphical clarity. 


Table 4.4. Wavelengths and Directions in Area IV 


Transform Area L(m) € (degrees) 
A27 970+65 179+0.5 
A30 845+49 1 7350. 5 
'836*48 182*0.5 
A31 7994-46 178+0,0 
784+44 REO. 5 


Depth ranges and effective median depths were 
predicted for area A30 using OFT data from areas A30 and 
A27, and for area A3l using OFT data from areas A30 and A3l. 
The depth range in area A27 is 150 to 185 m with an 
effective median depth of 170 m, and the range in A30 is 50 
to 190 m with an effective median depth of 135 m. Thus, 
using 970 m and 836 m for the measured wavelengths in areas 
A27 and A30, respectively, a predicted depth range of 102 to 
120 m with an effective median depth of 113 m was calculated 


for area A30. When compared to the actual effective median 
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Figure 4.7. 
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depth, the predietionehas sr l6=pemcenemenror. The observed 
change in wave direction is consistent with the bathymetry. 

The depth range in transform area A3l is 30 to 175 m 
with an effective median depth of 102 m. The predicted 
depth range is 44 to 151 m with a predicted effective median 
depth of 113 m, an error of ll-percent. The small error in 
the predicted depth range for area A31 is likely due to the 
large spatial overlap between the two transform areas. 

The waves in this area are propagating alongshore, 
rather than onshore and the wavelengths in transform area 
A31 may be due to the effects of shoaling which occurred 
somewhere to the north. Before shoaling occurred, off the 
image to the north, the deep water wavelength which 
refracted to that shown in A3l was presumably similar to 
that observed in the deeper water of transform area A27. 

5. Area V: Brightness Anomalies and Surface Gravity 

Wave Refraction 

Area V is located along the northern. edge of the 
dait, (Fig. 1.1). Figures 4.8, 4.9 and 4.12 show the 
imagery of the area from all three orbital passes. This 
area was covered by SAR images acquired during three 
different revolutions. As such, it is the only 
multi-temporal case in this study. One of the images (Fig. 
4.9) had such high intensity radar returns that visual 


interpretation was impossible. However, there are similar 
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Figure 4.8. Imagery from revolution 811 covering Area V, 
showing transform areas D17 and Dl6. 
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Figure 4.9. Imagery from revolution 289 covering Area V, 
showing transform areas Gl and G2. 
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Figure 4.10. Bathymetry of Area V, with transform areas Dl7 and DI6, 
(depths in meters, scale 1:250,000) 
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Figure 4.11. Bathymetry of Area V, with transform areas G! and G2, 
(depths in meters, scale 1:250,000). 
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Figure 4.12. imagery from revolution 323 covering Area V, 
showing transform areas Fl, F2 and F4. 
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Figure 4.13. Bathymetry of Area V, with transform areas Fl, F2 
and F4, (depths in meters, scale 1:250,000). 
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brightness anomalies in Figures 4.8 and 4.12 which appear to 
be related to the same bathymetric feature.. | 

There is an elongated brightness anomaly in Figure 
4.12 which corresponds to the shoaling area shown between 
58944'N and 58°45'N at 153°W longitude in Figure 4.10. 
Portions of a Similar-shaped anomaly can also be seen in 
Figure 4.9, although that anomaly is partially obscured by a 
large-scale cyclonic eddy feature. Winds are moderate and 
southeasterly; therefore, the mountains along the shoreline 
may be forcing part of the air flow along the shore, 
straining the surface radar cross-section such that a 
cyclonic brightness anomaly appears in the imagery. 

Figures 4.10, 4.11 and 4.13 show the corresponding 
bathymetry and transform areas and Table 4.5 lists the 
wavelengths and directions measured from the transforms. 
Note that transform D17 has two measured wave components, 
but Figure 4.10 shows only one for clarity. 

The relationship between measured wavelengths for 
transform areas Dl17 and Dl16 is reversed from the expected 
relationship. Figure 4.10 shows D17 is over shoaler water 
than D16, but has longer wavelengths. This anomaly is not 
explained by the uncertainty in the measurements; therefore, 
the transforms must be showing two separate wave structures. 


The spatial relationship of the two transform areas and the 
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shoreline makes it unreasonable to accept the two 


wavelengths as part of the same wavetrain. 


Table 4.5. Wavelengths and Directions in Area V 


Transform Area L(m) cx, (degrees) 

D17 902461 358+0.0 

891+58 353+0.5 
D16 588+24 346+0.5 
Gl 796+37 284+2,0 
G2 699+25 292+2.0 
Fl yov Y 4+0.0 
F2 1,229+106 19+1.0 
F4 1,106+88 E 5 


The predicted depth range for transform area G2, 
using OFT data from Gl and G2, is 19 to 49 m with a 
predicted effective median depth of 3l m. The actual depth 
range for Gl is 25 to 65 m with an effective median depth of 
40 m, and 30 to 75 m with an effective median depth of 35 m 
for G2, showing the predicted effective median depth has an 
ll-percent error. 

The change in wave direction from Gl to G2 is 
clockwise, MN ocn bathymetric contours oriented along a 
general northeast-southwest axis, which is consistent with 
the bathymetry in Figure 4.1l. 

The depth ranges in transform areas Fl, F2 and F4 
wea O Cto ds m, 50 to 125 m and 65 to 135 m, respectively 
with effective median depths of 32, 95 and 97 m, 
respectively. The predicted depth range for area Fl using 


data from F2 is 17 to 40 m with a predicted effective median 
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depth of 3l m, a 3-percent error. Data from transform area 
F2 were also used to predict a depth range and an effective 
median depth in area F4. The predicted depth range in F4 is 
40 to 99 m with an effective median depth of 76 m, an error 
of 22-percent. The differences in wavelengths and wave 
directions between transform areas Fl, F2 and F4 are 
consistent with the bathymetry, i.e., the wave directions 
are turning towards the normal to the bathymetry and the 
wavelengths are shortening as the wavetrain propagates 
shoreward. 
6. Area VI: Brightness Anomalies and Surface Gravity 

Area VI is located along the northwest shelf of 
Shelikof Strait (Fig. 1.1). Figure 4.15 shows the area's 
bathymetry and Figure 4.14 shows the imagery. Table 4.6 
lists the measured wavelengths and directions. | 

A small elongated patch of high reflectivity appears 
in the bay located in the upper right corner of the imagery 
(outlined area in Figure 4.14). This brightness anomaly 
CODEC directly to the broken, irregular bathymetry 
shown in Figure 4.15. There are also some brightness 
anomalies near the shoreline on the extreme left edge of the 
image in Figure 4.14 which correspond to the rough bottom 
shown in Figure 4.15. 

The dark areas which parallel the embayments to the 
north of the transform areas may be due to one of two 


conditions. These areas may be protected from the local 
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Figure 4.14. imagery from revolution 811 covering Area VI 
Eugsug transform areas C6; C7, C8. and C9. 
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Figure 4.15. Bathymetry of Area VI, (depths in meters, 
scale 1:250,000). 
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wind patterns (even though synoptic winds are southeasterly, 
making this an unlikely explanation); or they may be due to 
the influence of fresh water entering the strait from 


rivers. 


Table 4.6. Wavelengths and Directions in Area VI 


Transform Area Tom) SY (degrees) 
C6 838-44 326+2.0 
C7 658425 326+1.6 
cs 5:593 E] B2541.0 
C9 543-419 B23+1.0 


Depth ranges and effective median depths were 
predicted for transform areas C8 and C7. Data from area C9, 
with a depth range of 120 to 200 m and an effective median 
depth of 170 m, were used to calculate the predictions for 
C8. The predicted depth range for C8 is 108 to 153 m with a 
predicted effective median depth of 150 m compared to the 
actual depth range of 45 to 110 m and effective median depth 
of 80 m, an 88-percent error. 

The predicted depth range for:area C7 was calculated 
using data from area C6, which has a depth range of 140 to 
195 m and an effective median depth of 170 m. The predicted 
depth range of 75 to 99 m and predicted effective median 
depth of 85 m has a 6-percent error when compared to the 
actual depth range of 45 to 135 m with effective median 


depth of 90 m. 
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The measured wavelength differences between the four 
transform areas, C6, C7, C8 and C9, are consistent with the 
bathymetry shown in Figure 4.15. The very small changes in 
wave direction are consistent with waves propagating normal 
to the shoaling isobaths. 

7. Area VII: Brightness Anomalies and Density 

Stratification 

Area VII is located at the mouth of Kukak Bay (Fig. 
l.l). Figure 4.16 shows the imagery for the area and Figure 
4.17 shows the corresponding bathymetry. This area has the 
most dramatic brightness anomalies in this study. Because 
of the large modulations in the radar cross-section, this 
image is unsuitable for OFT analysis. Therefore, no attempt 
was made to interpret this image in terms of surface gravity 
wave refraction. 

There are a number of anomalies in the image in 
Figure 4.16 which appear to be related to bathymetric 
features shown in Figure 4.17. The largest anomaly 
(labelled A in Figure 4.17) is a sharp break in the image 
brightness, with high backscatter intensity on one side and 
very low intensity on the other. The bathymetry shows a 
small seamount less than 2.0 km upstream from the anomaly. 
The seamount sits at the mouth of a small submarine canyon 
with depths sloping from 250 to 150 m upstream of the 
feature and a fairly flat bottom with depths to 270 m on the 


downstream side. The feature has three peaks with charted: 
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Figure 4.16. Imagery from revolution 323 covering 
Area VII. 
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Figure 4.17. Bathymetry of Area VII, (depths in meters, 
scale 1:250,000). 
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depths of 170, 190 and 140 m respectively, progressing 
shoreward. 

These patterns could be the result of water column 
stratification due to fresh-water runoff. The cold, fresh 
Surface water layer would stabilize the atmospheric 
planetary boundary layer relative to the warmer, saltier 
surface water layer in the strait, producing a marked 
difference in the image brightness (due to the radar 
cross-section modulations) over the two different water 
masses. There are also internal wave signatures at the 
southeastern portion of the intensity break, the appearance 
of which is also characteristic of fresh-water runoff fronts 
where abrupt changes in stratification may occur (e.g., see 
Bice 34 of Fu and Holt, 1982). 

A second, smaller anomaly is located 500 m 
downstream from another, smaller seamount (marked B in 
Figure 4.17). The second seamount rises approximately 100 m 
from the base. The anomalous signature from feature B is 
offset from the Seamount, and resembles the pseudo-imaging 
of bottom relief similar to that found in Nantucket Shoals 
SAR imagery. 

Depths across Nantucket Shoals are generally less 
than 25 m, but the depths in this area range from 75 to 170 
m. This would indicate that the processes discussed in 


Section 11.B.3 may be valid for depths of 100 m or more, as 
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well as for the shallower depths found near Nantucket 
Shoals. 

There is a long, bright anomalous line located in 
the upper right-hand corner of the image in Figure 4.16. 
This anomaly (labelled C in Figure 4.17) is roughly located 
along a shelf break inshore of a narrow submarine canyon, 
and may be another example of a shelf break trapping a lens 
of cold, fresh-wateraunomu 

8. Area VIII: Surface Gravito Wave Ce aaan 
Area VIII lies outside Shelikof Strait, to the 


southeast of Kodiak Island and was chosen in an attempt to 


Table 4.7. Wavelengths and Directions in Area VIII 


Transform Area Em) S< (degrees) 
K3 866+38 229+2.0 
K4 923+60 179100 
K5 1,165+88 181+0.3 
K6 1,075+76 182+0.5 
K7 Ore 179 r05 
K19 886+55 176+0.5 
K20 Is Ta le FOS 
K21 1,388+131 164+2.5 


locate the channel seen in Figure 4.18 on the basis of wave 
refraction. A number of OFT's were attempted throughout the 
channel, and eight had distinct dominant wave numbers. The 
wavelengths and directions for the transform areas are 
listed in Table 4.7 and Figure 4.19 shows the imagery for 


the area. 
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Figure 4.18. Water depths in Area VIII, (depths in meters, 
scale 1:250,000). 
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Figure 4.19. Imagery from revolutson0c MN cV PS 
Area VIII, with transform locations. 


Figure 4.18 shows the wavelengths and directions in 
relation to charted Eee for the area. Transform areas K3 
and K4 have somewhat shorter wavelengths than the other 
areas and there is a general increase in wavelengths moving 
from left to right (from shallow water to deep) across the 
transform areas. However, the wavelength in K19 is 
equivalent to that in K3 and K4, and the water depth in Kl9 
is significantly deeper, according to the charted depths. 
The data in the diagram and Table 4.7 indicate a subsequent 
increase in depth moving from” left to right. 

Depth ranges and effective median depths were 
predicted for four transform areas, K6, K4, K20 and K19. The 
depth range in area K7 is 70 to 140 m with an effective 
median depth of 105 m. These values were used to predict a 
depth range of 50 to 96 m and an effective median depth of 
73 m in area K6. The actual depth range and effective 
median depth in K6 are 70 to 100 m and 85 m (14-percent 
error), respectively. 

The depth range and effective median depth for K4 
were predicted using the depth range in K5, 75 to 145 m with 
an effective median depth of 110 m. The predicted depth 
range in K4 is 46 to 84 m with an effective median depth of 
66 m compared to the actual range of 70 to 100 m with an 
actual effective median depth of 85 m (22-percent error). 

The depth range in area K21 is 100 to 160 m and the 


effective median depth is 130 m, which were used to 


103 


calculate the predictis range and the effective median 
depth in area K20. The predicted desc range form of 58 
to 88 m with an effective median depth of 73 m, compared to 
an actual 100 to 140 m range with an effective median depth 
of 125 m, has a 42-percent error. The charted data from 
area K20 was also used to calculate a predicted effective 
median depth and range for area K19. The K19 predicted 
values are 63 to 90 m with an.effective median depth of 77 m 
and the actual range is 100 to 150 meters with an effective 
median depth of 125 m, a 38-percent error. 

Since the depths are taken from a small-scale 
navigation chart (scale 1:2,100,000), it is possible that 
the K19 and K20 transform area wavelengths (shown in deep 
water) may indicate an uncharted extension of the shoal. 
Moreover, the change in wave direction between areas K2l, 
K20 and K19 would be more consistent if the 100 m depth 
curve were extended out into the channel at the locations of 


KI9vand su 


s SUMMARY 

Analysis of SAR images for eight case study areas in 
Shelikof Strait revealed many instances where significant 
bathymetric features were indicated by, 1) internal wave 
signatures, 2) radar cross-section anomalies, or 3) wave 
refraction detected from OFT's. Where possible, the 


mechanisms which caused anomalous patterns to appear in the 
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SAR imagery were described, and error budgets were 
calculated for the depths predicted from wave refraction. 
While there was one negative example, the analysis produced 
many positive examples. 

l. Internal Waves 

Three cases of internal wave signatures were 
observed in the SAR imagery of Shelikof Strait, and each 
signature was distinct from the other two. The internal 
wave propagating towards Augustine Island in Area I 
(identified by its SAR signature, a single wavetrain), is 
possibly generated by interaction of the internal tide with 
the shelf edge. 

A second internal wave signature appeared in the SAR 
imagery of Area II; a set of curved wave packets appearing 
to propagate outward from a single point source. The 
observed location of this internal wave signature is 
consistent with generation by interaction of the internal 
tide (during peak flood stage) with the seamount in the 
northeast entrance to Shelikof Strait. 

The third internal wave signature appeared in the 
SAR imagery for Area VII at the edge of a dramatic change in 
the backscatter intensity that may result from the presence 
of a fresh-water runoff lens. Density stratification occurs 
where there is a cold, fresh-water layer on the surface and 
the warmer, more saline ambient water lies on the bottom. 


When this occurs over a shelf with a steep slope into deeper 
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water, the self break tends to inhibit mixing of the layers, 
enforcing the stratification. Thus, interaction between 
water flow and the shelf break will tend to generate 
internal waves which appear near the edge of the fresh-water 
lens (Fig. 34 in Fu and Holt, 1982). Hence, the appearance 
of internal waves near the edge of a fresh-water lens would 
indicate a steep shelf break. 

25 Brightness Anomalies 

There were a number of cases where brightness 
anomalies in the SAR imagery could be related to either a 
shoaling bottom or a bathymetric feature. In Area II, a 
large anomaly appeared which, together with the internal 
wave signature discussed above, helped locate the large 
seamount near the northeast entrance of Shelikof Strait. In 
this case, the anomaly was backscatter return from surface 
waters TIN CCIV above the seamount. 

The brightness anomaly in Area IV corresponds to an 
area shown on the bathymetric map as having rocks awash. 
Thus, the increased brightness was likely due to modulation 
of the surface water's radar cross-section as well as to 
possible reflections from the rocks. 

In Area V, elongated brightness anomalies 
corresponding to a shoal area on the bathymetric map 
appeared in both SAR images of the area, and in Area VI 
there were several anomalies corresponding to a rough, 


irregular bottom topography. The increased backscatter 
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intensity from these areas is likely due to interaction 
between water flow and the rough bottom, modulating the 
Bragg waves and radar cross-section on the sea surface. 

Anomalies were expected to appear near Augustine 
Island (where rocks are known to break the surface), but did 
not (see analysis of Area I). There were also some cases 
where anomalies appeared, but the processes which produced 
them are not understood. 

In Area I, there are a number of rocky areas and 
rough bathymetry to the southeast of Augustine Island which 
should have modulated the surface radar cross-section and 
created brightness anomalies in the SAR imagery. However, 
no anomalies corresponding to these features appeared, 
possibly because the tide waters were too high for 
interactions between water flow and the bottom to modulate 
the surface radar cross-section. It is unclear, given the 
lack of information on prevailing oceanographic conditions, 
why this feature failed to produce an anomaly in the SAR 
imagery. This example clearly emphasizes that the absence 
of a brightness anomaly, or any other indicator of a shoal, 
in SAR imagery, can in no circumstances be interpreted as 
implying that an area is safe for navigation. 

Scale cre neS appeared in the imagery for 
Area III, which seemed to parallel the escarpment shown in 
the bathymetry. However, the lines are 2-10 km from the 


escarpment and the general water flow runs parallel to the 
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bathymetry. Thus, while the relationship between the 
ERR and the escarpment is not clear, the information 
can still be useful in locating bathymetric features. 

A similar type of bright line appears in the imagery 
for Area VII, situated roughly above a shelf break. The 
upper portion of the line appears to be the boundary of a 
fresh-water lens depicting the density stratification 
situation described previously. However, part of the lens' 
boundary crosses a canyon shown on the bathymetric map. The 
lower portion of the line does not appear to be associated 
with the fresh-water lens, but it is situated directly over 
the shelf break. Again, the mechanisms which produced the 
anomalous radar cross-section are not certain; however, the 
anomaly does appear to be related to the bathymetry and an 
interpretation of its presence can be useful. 

3. OFT Analyses 

Thirty areas in this study had a dominant wave 
spectra measurable by OFT (Table 4.8). Depth ranges and 
effective median depths were predicted for fourteen of the 
thirty areas. 

The predicted depth ranges and predicted effective 
median depths were calculated for a given transform pair 
using the OFT-measured wavelengths for each transform area 
and the depths shown in bathymetric maps of the area. The 
depth range and effective median depth for each transform 


area were determined by inspection of the bathymetry. In 
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Table 4.8 Optical Fourier Transform Data 
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Change in 
Transform Wavelength Wavelength Change in 
Pair Area (meters) (meters) Direction Direction 
B15 to I 1,825 18 
Bl4 378 1,447 348 30 
B15 to I 1,825 18 
B16 935 890 20 2 
A23 to II 1,160 175 
A22 845 315 180 5 
A26 to II 856 170 
A22 845 11 180 10 
A8 to ILII 2:002 189 
A4 1,150 1:272 181 8 
A5 to ITI 1,342 173 
A4 : 1,150 192 181 6 
A27 to IV 970 179 
A30 836 125 182 3 
A30 to IV 836 Py 
A31 784 52 177 4 
DI7 to V 902 358 
D16 588 314 346 po 
Gl to V 796 284 
G2 699 97 292 8 
F2 to V 1,229 ee 
Fl 722 507 4 T5 
F2 to V 1229 19 
F4 1,106 123 9 10 
C9 to VI 543 323 
C8 533 10 325 2 
C6 to VI 838 326 
C7 658 180 326 0 
K7 to VIII E O 179 
K6 1,073 196 182 3 
K5 to VIII 1,165 181 
K4 923 242 179 2 
K2l to VILE 1,388 164 
K20 1,068 320 167 3 
K20 to VIII 1,068 167 
K19 866 202 176 9 


the transform area having the longer wavelength (as 
determined from the OFT) the minimum and maximum depths were 
used in equation 4.2 to calculate estimates of the wave 
period. These wave periods were then used to estimate 
minimum, maximum and effective median depths associated with 
the other transform area. 

Transform areas Bl15/Bl4 and B15/Bl16 were excluded 
from the analysis because the large change in wavelength 
over such a short distance makes it improbable the 
OFT-measured spectra are from the same wavetrain. Also, 
transform areas D17/D16 were excluded because the areas' 
wavelengths were opposite from the tendency expected for 
onshore propagation of a wave group. In operational 
practice, this pair would be flagged as a probable 
"artifact" and set aside pending independent corroboration 
from subsequent SAR images. 

The predicted effective median depth was compared to 
the actual effective median range as determined from a 
bathymetric map and the error associated with the prediction 


was calculated by 


Error = 


median depth - median depth X 1003 


actual eff. predicted 
actual effective median depth 


The results of the predictions are given in Table 4.9. In 


general, the quantitative results are respectable. Ten of 
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Table 4.9 Optical Fourier Transform Prediction Accuracy 


Actual Predicted Predicted 
Actual Effective Effective Predicted Median 
Depth Median Predicted Median Depth Depth 
Transform Range Depth Shoaling Depth Range Error 
Pair (m) (m) Ratio (m) (m) q^ 
B15 to == <= == -~ -- =- 
B14 
B15 to => == == -- -— -- 
Bl6 
A23 to 225-— N59 145 
A22 110-160 1335 0. 53 70 62:76 48 
A26 to 100-160 PS0 
A22 110-160 135 0097 144 94-148 7 
A5 to 170-200 Io 
A4 50-200 150 0.73 130 SL 13 
A27 to P3O-—TO5 170 
A30 50-190 135 0.74 113 1025120 16 
A30 to 50-190 155 
A31 OS 102 0.88 1:13 44-151 14 
D17 to 
D16 a T == za ÉS - 
Gl to DOLO 40 
G2 SOS 35 1.00 31 19-48 11 
F2 to 50-125 95 
Fl 20-45 232 0235 21 17-40 3 
F2 to BUS 95 
F4 65-130 En oT 76 40-99 22 
CÌ to 120-200 170 
C8 45-110 80 0.96 Su 109—153 88 
C6 to 140-195 170 
C7 552135 90 0.62 85 15299 6 
K7 to 70-140 105 
K6 70-100 85 0772 73 50296 14 
K5 to 75-145 110 
K4 70-100 85 0.63 66 46-84 27 
K21 to 100-160 130 
K20 100-150 125 0009 13 58-88 42 
K20 to 100-150 12:5 
K19 100-150- i23 0.66 77 63-90 38 
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the fourteen predicted effective median depths had an error 
less than 30-percent, with a mean error of 24-percent. 
Predicted depth ranges, on the other hand, are 
consistently less than the actual depth ranges (Table 4.9). 
This occurs because the wave refraction estimated using 
OFT's selects a result consistent with the effective median 
depth for the area covered by each transform; the relatively 
small areas typically associated with shoalest and deepest 
depths in a pair of transform areas are simply not well 
represented by refraction of the dominant wave. Overall, 
the ratios of actual to predicted shoalest depths average 
1.16 and range from 0.42 to 1.77, and for deepest depths the 
ratios average 1.38 and range from 0.72 to 2.28. In cases 
where a reasonable prediction of median depth was obtained 
(30-percent or better), the ratios of actual to predicted 
shoalest depths averaged 1.07 and ranged from 0.42 to MEN 
and the ratios for deepest depths averaged 1.29 and ranged 


from La Tm So. 
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V. DISCUSSIONS AND CONCLUSIONS 


DES CUSS TON 

This study has explored the feasibility of using SAR 
imagery to locate hazards to navigation, thus enabling 
charting organizations to better utilize their limited 
resources. Analyses of anomalous patterns in SAR imagery of 
Shelikof Strait, Alaska, and wave refraction analysis using 
OFT-measured wave spectra were performed. 

Synthetic aperture radar imagery collections are 
generally not restricted by weather, cloud cover or hours of 
sunlight. Therefore, SAR sensors have an advantage over 
passive sensors, such as LANDSAT, for providing a 
quick-response capability to military charting 
organizations. As a military planning tool SAR can provide 
information on the safety of navigating through infrequently 
used waters, hostile areas, or waters which experience heavy 
cloud cover for much of the year, such as in South America. 
An analyst could provide information on these types of areas 
in a matter of days at a cost of hundreds of dollars, where 
it might take an advance survey party several weeks and tens 
of thousands of dollars to obtain the same information. 

Non-military charting organizations would also realize 
benefits by utilizing SAR imagery as a presurvey or planning 


tool. Where an organization has several areas to survey, 


eee 


and only enough funds to survey some of them, SAR images 
would reveal those areas with the most critical hazards to 
navigation, and thus requiring immediate attention. Again, 
SAR images would provide this information at a potential 
savings of tens of thousands of dollars. 

Anomalies occur in SAR imagery of ocean surfaces due to 
oceanographic and atmospheric processes modulating the sea 
surface radar cross-section. This study examined SAR 
imagery of Shelikof Strait, Alaska in an attempt to relate 
visual anomalies to bathymetric features and to locate shoal 
areas from OFT-measured wave spectra using wave refraction. 

A large seamount, a change in the ocean bottom slope and 
a shelf break were located using the signature in the 
surface radar cross-section left by the interaction between 
internal waves and these features.  Anomalous patches in the 
imagery were related to shoal areas, rocks awash, seamounts 
and freshwater lenses. Bright lines in the imagery were 
observed over or near escarpments and shelf breaks. 

In addition, linear wave theory was used to estimate 
shoaling ratios and depth ranges from OFT-measured wave 
spectra and wave refraction. The OFT-measured wavelengths 
and wave directions were obtained for locations in the study 
area most likely to experience measurable wave refraction. 
Predicted depths and shoaling ratios were compared to actual 
values determined from bathymetric maps, and an error 


calculated for the prediction. While the quality of the 


114 


results varied, in most cases the OFT-measured spectra 
qualitatively indicated the areas with the shoaler waters. 

This study also pointed out some risks involved with 
using the image analysis and the wave refraction analysis to 
locate hazards to navigation. In the image analysis portion 
of the study there were some areas where navigation hazards 
were known to exist, but which did not appear in the imagery 
(e.g., the rocky areas around Augustine Island in Area I). 
In the wave refraction analysis, there was one case 
(transform areas D17 and D16 in Area V) where the difference 
in wavelengths suggested the opposite sense of shoaling. 
However, this case can be flagged as a probable artifact 
without prior knowledge of the actual bathymetry on the 
grounds that it is unreasonable to expect wavelengths to 
increase during onshore propagation; in practice this OFT 
case would be flagged as an artifact pending confirmation by 
subsequent SAR images. Transform areas B15/Bl6 and B14/B15 
were also rejected on the apriori basis that it is unlikely 
a wavetrain could experience wavelength shortening of over l 
km over a distance of 5 km, without breaking. 

Caution and prudence must be exercised when interpreting 
the results of SAR analysis for detecting hazards to 
navigation. If an anomaly appears in the imagery, a 
bathymetric feature is likely to exist in the area, but the 
absence of a signature does not necessarily imply that an 


area is hazard-free. In SAR OFT wave refraction analysis, 
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observed changes in wavelength are not always correct 
indicators of shoaling; an analyst must also consider the 
magnitude of the change and other circumstances to ascertain 
whether the two spectra are likely to represent refraction 
of a single wavetrain. 

Independent knowledge of environmental conditions, such 
as tides, tidal currents and wind patterns is also necessary 
for correct interpretations of SAR images. Proper ground 
truth information, such as tidal heights, actual times for 
the tides, weather conditions and both surface and 
subsurface current data, would have significantly enhanced 
this study. The lack of concurrent oceanographic and 
meteorological data has imposed the need for speculative 
interpretations of many of the phenomena observed in the 
study. Future research projects and operational 


applications of this technique should include provisions for 


concurrent acquisition of these data whenever possible. 


B.  CONCCUSTONS 

The results of this study show that SAR images contain 
potentially useful information about shoals and other 
bathymetric features. An experienced, trained analyst could 
use SAR imagery to derive valuable presurvey information for 
planning hydrographic surveys. He could identify potential 
hazards to navigation, determine the general shape of the 


bathymetry using SAR signatures, and identify areas and 
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approximate depths of shoaling waters using OFT analysis to 
measure wave refraction. This information could then be 
used to assign priorities to areas to be surveyed. 

An analyst should not, however, expect to assign 
accurate depths to any areas based on either visual 
information or OFT analysis alone. Furthermore, while SAR 
images often indicate navigational hazards, the absence of 
anomalies does not mean there are no hazards present. 

This study also shows that wave refraction can be 
measured by OFT's well enough to warrant further study, 
albeit in an area where the refraction can be more easily 
measured. Most of the areas studied here are sheltered from 
the Pacific swell, and the waves illuminating the 
bathymetric features were locally wind generated waves of 
presumably smaller amplitude and wavelength. Shelikof 
Strait is, therefore, perhaps not the ideal test area for 
exploring the application of SAR-measured wave refraction 
(using OFT's) to detect shoals; a study of Pacific atolls, 
which are illuminated by relatively large amplitude swell 
generated by storm winds with virtually unlimited fetch is 
likely to be more conclusive. From this standpoint it would 
have been better to have evaluated SAR imagery from such 
reasonably well-charted Pacific islands as Guam, Yap, 
Saipan, Iwo Jima, Okinawa, and various islands and 
reef/shoals in the Hawiian Archipelago. These shoals are 


well illuminated with Pacific swell, and many areas are 
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reasonably well surveyed. In addition, such a study would 
have been valuable in highlighting possible uncharted shoals 
in subareas which have not been recently surveyed. 

Despite the difficulties in analyzing a sheltered area 
such as Shelikof Strait, there is clearly valuable 
information to be gained through SAR imagery analysis. A 
trained analyst who is experienced in interpreting 
oceanographic and meteorological phenomena, and who also 
understands the radar imagery, could develop a presurvey 
bathymetric feature analysis which could lead to significant 
cost savings in survey implementation. In addition, this 
analysis would enhance the safety of survey operations by 
alerting shipboard personnel to locations of possible 


navigational hazards. 


C. RECOMMENDATIONS 

The results of this study lead to several 
recommendations for further work in SAR imagery analysis and 
for utilizing that analysis. Additional research is 
necessary to better determine the accuracy of depth 
inferences from wave refraction measured by OFT's of SAR 
imagery. This research should be in an area open to the 
ocean, where the OFT can track wavetrains as they change 
from large amplitude, deep-water swell through transition 
waves to shallow-water waves. Also, the transform areas 


should overlap, such as transform areas A30 and A3l in Area 
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IV, rather than having discrete transform areas, as was the 
case throughout most of this study. This would provide a 
better tracking capability and possibly provide information 
on exactly where different wavetrains begin to mix such that 
there is no dominant wave spectra. 

A second recommendation is that further work be done to 
evaluate systematic procedures for subjective and 
Were quantitative interpretation of brightness anomalies and 
internal wave signatures. These methods should also 
incorporate information on local oceanographic, 
WeEeorologrcal, and hydrological (i.e. river runoff) 
conditions to allow less ambiguous interpretations of 
certain features, and to allow a better discrimination of 
possible wind stress anomalies than was possible in this 
study. 

l The third recommendation is that existing SAR data be 
utilized in an interpretive mode as part of the operational 
presurvey planning process. Image analysts with experience 
in both oceanographic and radar phenomena can extract 
valuable information from SAR imagery, and both interpretive 
methods and procedures for using such information in 
planning surveys can be best developed in an operational 
environment. Once surveying and charting organizations make 
a committment to utilizing SAR imagery, the methods and 
procedures developed should lead to significant improvement 


in efficiency and cost-effectiveness of survey operations. 
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APPENDIX A: OPTICAL FOURIER TRANSFORMS 


A. APPLICATION TO WAVE REFRACTION 

Experiments by Abbe and Porter (Goodman, 1968), showed 
that coherent, collimated light passing through an object 
transparency and then through a lens separates the object 
pattern into its Fourier components at the transform lens' 
focal plane (Fig. A.1). A simple illustration of the 
transform is shown in Figure A.2. If the transparency has a 
pattern similar to the one in Figure A.2a, then the Fourier 
components shown in Figure A.2b will appear at the focal 
plane. The distance from the origin of the axes to the 
first dot is the k-space wavelength, or wave number, of the 
fundamental harmonics. .It also follows that the distance 
from the origin to each succeeding dot is representative of 
the next higher harmonics. If the pattern shown in Figure 
A.2a is rotated as shown in Figure A.3a, the subsequent 


transform also rotates as shown in Figure A.3b. 


B. GENERATING THE OFF 
Optical Fourier Transforms (OFT), were generated using a 
helium-neon laser mounted on an optical bench (Fig. A.4). 


The laser beam is expanded using a pinhole device and 
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Figure A.2. OFT of evenly-spaced grating. 
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Figure A.3. OFT of evenly-spaced grating after rotation. 
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collimated using a lens placed so the pinhole is at its 
focal pointa 

The distance between the two lenses must allow a roll of 
film and a plate containing the focusing aperture to pass 
between them. The transform of an object on a film 
transparency occurs at the focal plane of the second lens 
(Fig. A.l). To focus the laser beam on a particular area of 
the imagery, a metal plate witha 1 S aperture was 
inserted between the collimating lens and the film imagery. 
With the scale of the imagery at 1:500,000 it was possible 
to determine the dominant wavelength and wave direction for 
any given 25 "n area. 

The transforms were displayed on a piece of white 
cardboard, (with calibration marks), placed in the focal 
plane and then photographed. A small hole in the cardboard 
was aligned with the center of the transform to allow the 
intense light of the center beam, and very long wave 
components, to pass through; if this was not done, the 
desired transforms could not be distinguished from 


background brightness in the photographs. 
C. TRANSFORM MEASUREMENTS 


The transform wave number components were measured using 


a zoom transfer scope. Using the geometry of Figure A.95, 
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Figure A.6. 
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where dx is the measured image distance between 2p and the 


origin, dy is the measured image distance between F1 and Fa. 


The scaling factors l/dx and 2/dy corrected the coordinate 
measurements relative to a true dx of l cm and a true dy of 


2 cm. anus; 
s oe 

and the wavenumber in k-space, 
E MI 2r; 

and the transform heading, 


a= ta eon 


From Figure A.6, the surface wave Structure heading, 


where a' is the satellite heading and the sign is determined 


by the quadrant of the transform heading. In determining 
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a", it was assumed that surface waves were always 
propagating shoreward into shoaler water, (with the 
exception of Area II, section IV.B.2, where the surface 
waves were assumed to be propagating in a southerly 
direction. 

Four optical fibers were attached by one end to a red 
light source and the other end placed in the cardboard to 
produce fiducial marks in the transforms. The ends of the 
fibers were placed such that each dot was one centimeter 
from the center hole, and such that if two lines were drawn 
between each pair of opposing dots, the lines would be 
orthogonal, creating a rectangular coordinate system with 
the small hole as the origin. The fiducial marks served to 
correct for geometric distortion caused by the camera's 
angle relative to the focal plane, and for other systematic 
distortions as well. 

Photographs were taken using a GAF model L-14 35mm 
camera with extension tubes between the aperture and lens. 
The extension tubes allowed the transform display to be 
mapped onto the camera film at a ratio of almost one-to-one. 
Otherwise, it would have been necessary to enlarge each 
photograph in order to make measurements, which would have 
introduced another error source. 

Kodak HIE (High Speed Infrared Film 2481) was used in an 
entirely darkened room to photograph the transform produced 


by the red-light laser. The aperture and exposure time 
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varied according to the brightness of the transtorn wa eles 
the photographs were developed, a zoom transfer scope was 
used to enlarge them slightly for each measurement. 
Coordinates in an x,y rectangular coordinate system were 
then measured for each pair of transform component dots. 

The distance from the center hole to each fiducial mark 
was also measured to correct for distortions caused by 
camera angle, transform lens and the zoom transfer scope. 
The scaling factor to correct that measurement to a true l 
cm on both x- and y-axes was applied to the measured 
coordinates of the transform component dots. 

Wavenumber measurement errors were oni the order NE ne 
mm. This error was a function of the repeatability of each 
measurement, and is therefore, independent of the actual 
measured distance. For a measured k-space distance of 1 mm, 
there is a probable error of 10-percent which, from equation 
A.2, converts to a surface wavelength error of +188 m (L = 
1880 m). For a k-space measurement of 20 mm, there is a 
probable error of 0.5 percent, which yields a surface 
wavelength error of +0.5 m (L = 94 m). Thus, the shorter 
the measured k=space distance these eae: pet Mempe Weemiad@c mcm 
measurement error. 

There are other error sources present on the transforms 
in the form of distortions due to the quality of the lenses 
and the lack of precision in the focusing aperture that were 


used to generate the OFT. Therefore, a direct mapping of 


2 


wavelength L on the imagery to the transform is not precise. 


Since the wavelength on the imagery can be expressed as, 


A.l L = 2m/k, 


where k 1S the wavenumber, we may express the wavelength on 


the transform as, 


A.2 L = K/(sr), 


where K is oceanwave image scale factor (m/mm), s is a 
dimensionless transform coefficient and r is the measured 
distance of the transform component from the center hole 
(mm). Note that the coordinate system was changed from 
rectangular to polar coordinates for ease of computation in 
determining wave direction. A transform for a grating of 
known wavelength was generated, and the distance to the 
fundamental transform component was measured. For K=l, the 
known values of L and r were substituted in equation A.2 to 
obtain the value of 0.266 for s. However, for the scale of 


the imagery at 1:500,000, the value of K is assigned as 500. 
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APPENDIX B: UOSURVEPDEHBESIBORDES 


This appendix summarizes dates, scales and line spacing 
information on hydrographic surveys in the Shelikof Strait 
study area. 

Figure B.l shows the geographic locations of the 
bathymetric maps for Iliamna, Mount Katmai and Afognak used 
in this Seud: 

Figure B.2, in conjunction with Table Bl gives the 
survey history for the bathymetric map covering the Afognak 
area. The labelled area in Figure B.2 is cross-referenced 
with the appropriate survey date, scale and line spacing to 
determine the quality of the bathymetric map. 

In a similar manner Figure B.3 is cross-referenced with 
Table B.2 and Figure B.4 with Table B.3 to determine the 
survey quality of the Mount Katmai and Iliamna areas, 


respectively. 
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